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Abstract 

The transverse velocity of a spiral galaxy can be measured to an accuracy 
~ 60 km s~"^ by making parallax observations of quasars being microlensed by stars 
in the disk of the galaxy. To make the measurement, a quasar must be located 
behind the disk of the galaxy between about 1 and 2 scale lengths from the center. 
The quasar must then be monitored for microlensing events and the events followed 
simultaneously from the Earth and a satellite in solar orbit, preferably at ~ 30 AU. 
A systematic search in a volume within 21, 000 km s~"^ could locate quasars with 
i? < 23 behind a total of ~ 1900 galaxies. The rate of lensing events (and hence 
galaxy velocity measurements) would be ~ 3yr~"^. The events would have typical 
characteristic times lo~^ ~ 3yr. Under the assumption that the mass spectrum of 
lensing objects is the same in other spiral galaxies as in our own, the observations 
could be used to measure the Hubble parameter to an accuracy of 5%. 

Subject Headings: astrometry - gravitational lensing 
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1. Introduction 



Measurement of the transverse velocities of galaxies could provide important 
information about the large scale structure of the universe. For example, Dekel, 
Bertschinger, & Faber (1990) have shown that one may reconstruct the local galaxy 
velocity field (including the transverse components) from measurements of peculiar 
radial velocities, provided that one assumes that the velocity field is irrotational. 
By measuring the transverse velocities one could test this assumption. Several 
workers have attempted to measure the bulk fiow of the nearby universe out to 
various redshift limits from 3000 to 15,000 kms~"^ (Courteau et al. 1993; Mould 
et al. 1993; Lauer & Postman 1994). In particular, they have sought but so far 
failed to find a volume sufficiently large that it is stationary with respect to the 
cosmological reference frame defined by the microwave background radiation. The 
measurement of transverse velocities would yield an independent determination of 
these bulk fiow velocities. 

In addition, all peculiar radial velocity measurements are based on the com- 
parison of a redshift with a distance as determined from a distance indicator. This 
approach has two intrinsic drawbacks. First, since the distance errors are pro- 
portional to the distance, the fractional errors in the peculiar velocity are also 
proportional to the distance. Thus, even if the distance errors are only ~ 6% as 
claimed for the method of surface brightness fiuctuations (Tonry 1991), the er- 
rors rise to order unity for 600kms~"^ peculiar velocities measured at redshifts of 
10, 000 km s~^. Other methods, such as TuUy-Fisher for spirals and Dn-cr for ellip- 
ticals have substantially larger errors and hence smaller effective ranges. One would 
like to have a method of measuring velocities for which the errors are independent 
of distance. Second, all secondary distance-indicator measurements are based on 
correlations between various observed properties of galaxies. The physical bases 
of these correlations are incompletely understood. One must therefore worry that 
subtle selection effects might introduce systematic errors in the distance determi- 
nations. Fven small systematic errors in the fractional distance measurement can 
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be serious at large distances where the pecuhar velocity is itself a small fraction 
of the redshift. One would like to have a method of measuring velocities that is 
based directly on known physical principles. 

Here I present a method of measuring the transverse velocities of galaxies from 
the parallax of microlensing events generated by stars within the galaxy. The 
method is based on simple and well-established physical principles. For spiral 
galaxies, the accuracy of the measurement is limited only by the dispersion of stars 
in the disk, typically ~ 40kms~"^ in each transverse direction. In particular, the 
limiting errors are independent of distance. The measurements can be carried out 
using existing technologies. They do however require a considerable investment of 
resources including a ~ Im telescope in solar orbit, preferably at ~ 30 AU, as well 
as an intensive search for faint quasars behind thousands of spiral galaxies. 

Grieger, Kayser, & Refsdal (1986) were the first to suggest that the transverse 
velocities of distant galaxies could be estimated using gravitational microlensing 
of distant quasars. They pointed out that if a quasar is highly magnified by a 
lensing galaxy, then the stars in the galaxy produce a network of microlensing 
caustics. If the line of sight from the Earth to the quasar should pass over one of 
these caustics, the magnification would reach very high levels. By comparing the 
times of maximum magnification as seen from the Earth and a satellite, one could 
measure the "reduced" transverse speed up to an unknown geometric factor. The 
reduced velocity v is related to the true transverse velocity v of the lens relative 
to the observer-source line of sight by 



where -Dqs, -Dls, and -Dql are the angular diameter distances between the observer, 
source, and lens. 

I have previously shown that the reduced transverse velocity of Massive Com- 
pact Objects (MACHOs) can be measured by comparing the microlensing light 
curves as seen from the Earth and a satellite (Gould 1992, 1994a, 1994b). Note 
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that for MACHOs both components of v can be measured without ambiguity, 
whereas for the highly magnified quasar only the reduced speed can be measured 
and then only up to an unknown correction. This important difference is rooted 
in the topologies of the caustic structure of each type of lens. MACHOs are iso- 
lated point masses and so induce a circularly symmetric magnification field about 
an isolated point caustic. The difference in light curves therefore gives two di- 
mensional information about the motion of the lens relative to the Earth-satellite 
separation vector (see Figs, f and 2 in Gould f994a). By contrast, the highly mag- 
nified quasar field is laced with linear caustics. Hence, information about relative 
motions is inherently one dimensional (see Fig. f in Grieger et al. f985). 

The method I propose here is simply to search for non-highly magnified quasars 
behind the disks of spiral galaxies. The stars in such galaxies give rise to isolated, 
circularly symmetric magnification fields. The method developed for measuring 
MACHO reduced velocities therefore applies directly to such lensing events. In the 
case of distant quasars lensed by stars in relatively nearby galaxies, (I-|-2;)_Dos /-Dls 
is close to unity and can in any event be accurately measured. Hence the transverse 
velocity can be inferred from the measurement of the reduced velocity. This situ- 
ation is much superior to the MACHO case where the ratio Dqq/Dj^q is generally 
unknown. 

The reader who is interested in obtaining a clear picture of the procedure for 
determining reduced velocities from parallax measurements is urged to study Gould 
(1994a, 1994b). 

I propose an experiment that would monitor quasars brighter than i? ~ 23 
behind ~ 1900 galaxies over ~ Stt steradians out to a redshift ~ 21, 000 km s~"^. 
I estimate that one could measure ~ 3 transverse velocities per year i.e., ~ 100 
over three decades. Such a large sample of high precision measurements would 
both serve as an independent check on models derived from radial measurements 
and greatly enhance the accuracy of the models. For example, assuming that the 
galaxies themselves have random motions ^ 200kms~"^ in each direction, such 
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a sample would measure the bulk flow of the sampled volume with a precision 
^ 25kms~^. The bulk flow of each of fO subregions could be measured to ^ 
80kms~"^. These determinations would be free from the "rocket effect" (Kaiser 
f987) that afflicts studies based on radial peculiar velocities. 

In addition, the experiment would have several other important applications. It 
could be used to measure the Hubble parameter, to determine the ratio of stellar 
mass to light of spiral disks, and to study the structure of quasars on scales of 
100-1000 AU. 



2. Lensing Events 

Suppose that a distant quasar is lensed by a star in a relatively nearby galaxy 
at a distance D = -Dql- I assume that the optical depth r for lensing by stars 
along the line of sight through the galaxy is low 

47rG'Si)oLi)LS MJW 

where E is the column density of stars and where I have assumed Dj^q/Dqq ~ I. 
Equation (2.1) can be rewritten as vrE/Ecrit ^ 1 where Edit is the critical density. 
Hence r ^ I is equivalent to the surface density being well below the critical level 
required for macrolensing. In turn, this usually implies that the shear is low. In 
any event, I will assume that this is the case. Hence, the star can be treated as 
an isolated microlens and the formalism that I have previously developed (Gould 
1994a) can be applied directly. The Einstein ring radius, rg is given by 



re = — = 280 AU — , 2.2 



where M is the stellar mass and -Dcoma is the distance to the Coma cluster which I 
take to be 100 Mpc (i.e. Hq = 70kms~"^ Mpc~"^). The Einstein ring projected onto 
the source plane has size {Dos/DoL)re ~ 10^ AU (M/O.IM0)l/2(i)/i)coma)"^/^ 
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where I have assumed that the quasar is at a cosmological distance. This size is 
much larger than the size of the quasar continuum source. Therefore, quasars that 
are much brighter than their hosts can be treated as point sources. The lensing 
event as viewed from the Earth can be described by a three parameter hght curve 

AW)] = ^p^^; <t) = ^uHt-t,y + (3^ (2.3) 

where A[x) is the magnification, x is the projected separation of the lens and 
source in units of the Einstein radius, lo~^ is the characteristic event time, to is 
the time of maximum magnification, and /3 is the dimensionless impact parameter. 
The light curve as seen from a satellite at position r relative to the Earth has a 
similar form but with different parameters /3' , t'^ and cu' ~ cu. One can form a 
vector Ax = (cuAt, A/3) where |A/3| = it /3| and At = t'g — to. Ignoring for the 
moment the four-fold ambiguity in A/3 (see § 5.3), one can determine the transverse 
velocity by measuring the parameters of each light curve and using the relation 
V = c<;rAx/(Aa;)^. Note that Ax = r/r-g should be in the range O.I ^ Ax ^ I 
in order to make a reasonably accurate (^ 10%) velocity measurement. Hence, 
to measure the transverse velocities of galaxies out to three times the distance of 
Coma, the satellite orbit should be at least ~ 30 AU. The characteristic time of a 
typical event would be 

^ ^ ^ 2.2yrf^V7^V7 / V'. (2.4) 
V \0.lMcJ \DcomJ V600kms-V ^ ^ 
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3. Proposed Experiment 

In order to measure the transverse velocity of a spiral galaxy, one must first 
find a quasar behind the galaxy's disk. The number density of quasars brighter 
than Bo = 22 is ^22 ~ fOOdeg"^ (Hartwick & Schade f990). The number brighter 
than Bo = 23 is ^23 ~ f50deg~^ (P. Osmer f994, private communication). To 
estimate the number of galaxies with quasars behind them, I will assume that all 
galaxies are face on and that quasars with Bo < 23 can be detected between I and 
2 disk scale lengths. I take M3I to be a 'typical' L* spiral, and will assume that 
there are 1.2 X IO~^/i^Mpc~^ ~ 4 X IO~^Mpc~^ such galaxies (Binney & Tremaine 
1987) with scale length h = 6.4 kpc and total blue luminosity oi L = 2.4 X 10"^*^ Lq 
(van der Kruit 1989). I take the stellar mass to blue light ratio to be A Mq / Lq^ 
implying a surface stellar density at radius £ of 

E = 375exp(-£//i)M0pc-^ (3.1) 



The adopted assumptions are obviously rather crude. I justify this approach below. 

I then find that of galaxies at a distance D a fraction / = 3Tr[h/D)^n23 = 
0.0I9(-D/-Dcoma)~^ have quasars behind them between I and 2 scale lengths. For 
those galaxies with quasars at radius r = 0.023 exp(— £//i)(_D/_Dcoma)- Hence 
the mean optical depth over the disk is r = 0.005 (-D/Dcoma)- The rate of lensing 
events by an average galaxy with a quasar behind it is then F = {2Itt)tlo or 

F = 1.4 X 10"^ yr"^ . 3.2 

^ \0.1MqJ \Dcou.J V600kms-V ^ ' 

Consider now a volume limited sample restricted to the Sir steradians [\b\ ^ 
20°) within the galactic polar caps and inside a distance -Dmax- The total number 
of L* galaxies in the sample is = I2500(-D/-Dcoma)^- Of these, a total of 

iVQ = 700 ^^, (3.3) 

i^'Coma 
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will have quasars behind them. The rate of lensing events will then be 

In deriving these results I assumed that all galaxies are face on. If the galaxies 
are at all inclinations but all the quasars behind the galaxies are located, the 
total number of quasars will be reduced but the optical depth for each will be 
proportionately increased. Hence, equation (3.4) remains valid. Of course, if a 
galaxy is highly inclined, the extinction through the inner disk will be increased 
and some of the background quasars will be lost. On the other hand, the disk will 
then have a relatively high optical depth even at several scale lengths so that it 
would be profitable to search a greater part of the disk for quasars. I also assumed 
that Bq = 23 quasars could be found even at 1 scale length where the extinction 
might be ~ 1, so that B ~ 24. As I will discuss in the next section, finding and 
monitoring quasars that faint may prove difficult. On the other hand, the number 
of quasars Bq < 22 within the inner scale length will not be insignificant and these 
(with their higher optical depth) will compensate for the lost faint quasars further 
out. At low latitudes some of the fainter quasars will be lost due to foreground 
extinction, but this results in only a ~ 10% reduction in the sample averaged over 
all angles. Because the quasar luminosity function is not rising rapidly at these 
apparent magnitudes, the estimates that I have made are not very sensitive to the 
exact cutoff assumed. 

From equations (3.3) and (3.4) and taking account of the 10% reduction due to 
foreground extinction, I find that a survey out to -Dmax = 3-Dcoma (~ 21, 000 km s~"^) 
would monitor Nq ~ 1900 quasars and measure the transverse velocities of ~ 3 
galaxies per year. After 30 years, ~ 100 velocities would be measured, implying 
a mean sampling scale of ~ 6600 km s~"^ over the whole volume. The measure- 
ment errors would be (9(10%) in each direction (see § 5.3). For typical expected 
transverse velocities ~ 400kms~"^ in each direction, this error is comparable to the 
intrinsic error arising from the rms difference between the lensing star's velocity 
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and the mean velocity of the spiral disk at its location. Hence the overall errors 
would be ~ 60kms~"^. This is less than the expected random motions of galaxies 
relative to the mean peculiar velocity flow. 

4. Scientific Byproducts 

Here I briefly describe the scientific questions which could be addressed with 
such a survey other than the results pertaining to large-scale structure described 
in§f. 

4.1. Measurement of the Hubble Parameter 

Each transverse- velocity measurement will also yield rg, the radius of the Ein- 
stein ring induced by the lensing star. Erom equation (2.2) and Hubble relation 
one finds M = r^c/^GHoz where Hq is the Hubble parameter and z is the galaxy 
redshift. Hence the mass spectrum of the lensing stars in an ensemble of galax- 
ies can be measured up to an overall factor of Hq. The errors in the individual 
mass measurements would have errors (~ 20%) that are small compared to the 
half- width of the mass spectrum (which I take to be ~ 50%). The mass spectrum 
of lensing stars in the Galactic disk can be determined by a similar experiment 
(in many ways a small prototype for the one proposed here) that would measure 
parallaxes of MACHO events seen in ongoing lensing experiments (Gould 1994a, 
1994b, Han & Gould 1994). Hence, if one assumed that the Galactic spectrum was 
the same as the composite extragalactic spectrum, the Hubble parameter could 
be measured to ~ 5% from 100 events. Of course, the mass spectrum may vary 
according to galaxy type or position in the galaxy. However, the variation of the 
spectrum as a function of galactocentric radius in the Milky Way can be measured 
from MACHO parallaxes (Han & Gould 1994), and a similar test can be performed 
by dividing the extragalactic sample into subsamples according to the radial posi- 
tion of the quasar line of sight. Similarly, the extragalactic sample can be divided 
up into subsamples of different galaxy types to test for different mass functions. 
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4.2. Ratio of Stellar Mass to Light 

The experiment allows one to measure the mean value of the surface mass 
density along all the lines of sight to program quasars by combining the observed 
mean optical depth with the galaxy distance distribution. This can be compared 
directly with the mean surface brightness to get a ratio of stellar mass to observed 
light. One could also find the correction for self absorption by studying the optical 
depth versus surface brightness as a function of galaxy inclination. 

4.3. Central Regions Of Quasars 

As I discussed in § 2, the Einstein ring projected onto the source plane will have 
a radius of order ~ 10^ AU, so that the form of the lensing event will generally not 
be affected by the much smaller structure of the quasar source region. However, 
a fraction of events will have impact parameters fi < fi^ and hence will come 
within f 0^/3^ AU of the central source. Thus over 30 years, f events will come 
within fOOO AU and f within fOO AU. These regions can then be probed with 
ground-based spectra and or analysis of the deviation of the light curves from the 
expected microlensing form. 

5. Practical Requirements 

Here I give a brief description of some of the potential obstacles to carrying 
out the experiment and how these obstacles might be overcome. 

5.1. Quasar Identification 

The surface brightnesses of spiral disks are fainter than the sky, so it is no more 
difficult to locate a quasar behind a galactic disk than it is to find an object of 
the same apparent magnitude in the field. However, there are potential problems 
posed by contamination. 
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Possible contaminants are foreground stars, faint galaxies, and stars, globular 
clusters, and HII regions in the spiral galaxy. Quasars with 2; < 2.2 generally have 
a substantial UV excess. If this is the principal method by which the quasars are 
selected, then virtually all globular clusters, HII regions, and stars in the galaxy 
will be excluded. All Galactic foreground stars except young disk white dwarfs and 
distant horizontal branch stars will likewise fail the test. Horizontal branch stars 
at i? ~ 23 are so distant (~ 250 kpc) that they are extremely rare. Extremely hot 
white dwarfs are also rare and moreover, can be found by their parallax offset as 
seen from the Earth and the satellite. Most background galaxies can be eliminated 
by the color criterion and most of the rest can be resolved in good ground-based 
seeing. Experience with the Hubble Space Telescope (e.g. Bahcall et al. 1994) shows 
that a few hours of observations on a Im satellite telescope could easily resolve the 
remaining ambiguities. These observations could be undertaken as the satellite 
traveled toward its Neptune like orbit. 

In any event, contamination at the level of a few tens of percent is not critical. 
Ground based spectra could be taken of the brighter sources and the fainter ones 
could simply be left in the sample and monitored over the first few years for 
variation. 

5.2. Event Identification 

Some indication that an event may be going on will come from the light 
curve (which to first order will be the same as seen from the Earth and the 
satellite). However, since quasars are intrinsically variable, this signature would 
not be definitive. Absolute proof of lensing can come only by demonstrating 
that the ground-based and space-based light curves are different. Regardless of 
any intrinsic variation, the fractional difference between the two curves will be 
AA = 8{r /re) cos 9 /[x^x"^ -\- 2){x'^ -\- 4)] where 9 is the angle between the pro- 
jected Earth-satellite and source-lens vectors. Since r/rg ~ (9(10%), this effect is 
AA ~ (9(2%) at x ~ 1.3. Since events have characteristic times lo~^ ~ 3yr, the 
effect would be noticed provided that ~ 3% photometry were carried out from the 
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satellite once every other month. Once it was noticed, a definitive set of three 1% 
measurement could be made to confirm the lensing character of the event. 

5.3. Measurement of Velocity 

An important difference between measuring the reduced velocities of a Galactic 
MACHO against a background star and of an extragalactic lens against a back- 
ground quasar is that the latter source is variable. Hence, the overall light curves 
as seen from both the Earth and the satellite deviate from standard microlensing 
light curves and the accuracy with which the basic parameters cu, /3, and to can 
be measured from these curves may be significantly compromised. However, the 
ratio of the fiuxes seen by the Earth and satellite is independent of the quasar's 
variability. Moreover, it is possible in principle to reconstruct cu. At, /3 and /3' from 
the ratio of light curves. Unfortunately, for a given level of photometry, the errors 
are then much greater than if one could measure the unadulterated microlensing 
curves directly. 

To investigate this problem quantitatively, I perform the following calculations. 
I assume that events are typically recognized at separations of f .3 Einstein radii. 
Before this time, I assume that the satellite makes 3% photometry measurements 
every other month and afterward it makes 1% measurements every week (as re- 
quired to break the degeneracy among the 4 possible solutions - see below). The 
ground-based measurements are assumed to have substantially smaller errors. I 
consider the relatively difficult case of a galaxy at the edge of the volume, with 
r-g = 350 AU and lo~^ = 4yr. I then estimate the errors in measuring cu, to At, 
and A/3-1- under two sets of assumptions: first, that some information about io and 
/3 are available from the overall light curve and second, that the parameters are 
estimated only from the ratio of the two curves. I find that for /3 ^ 0.5 and for 
many orientations of v relative to r, it is possible just from the light-curve ratio 
to measure Ax to ^ fO% and io to ^ 20%. The latter value then dominates the 
errors. However, the situation deteriorates sharply if v • r ~ 0. In addition, for 
/3 ~ I the errors rise to order unity at all orientations. On the other hand, if the 
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distortions of the overall light curve due to quasar variability are not severe so that 
the time scale lo~^ can be estimated to ~ 10% then the typical errors even at /3 ~ 1 
are only ~ 15%. 

At present therefore, the largest single uncertainty in outlining the proposed 
experiment is our poor knowledge of the statistics of quasar variability. After 
the MACHO project (Alcock et al. 1993, K. Cooke 1994, private communication) 
identifies the quasars in their Large Magellanic Cloud fields, it will be possible 
to make a systematic study of quasar variation and to estimate the effects on 
the galaxy velocity measurement. Here I will give only a brief overview of the 
velocity determination. I will assume (as I believe to be the case), that quasar 
variability does not qualitatively compromise the measurement of the overall light- 
curve parameters. 

The accuracy of the velocity measurement then depends on the accuracy of 
Ax = (c<;At,A/3). The difficulty of determining Ax is greatly reduced by the 
fact that typically r << rg. This means that for D ^ -Dcoma both light curve 
trajectories usually pass on the same side of the lens. That is, |A/3| = A/3_ (see 
Fig. 1 from Gould 1994a). For example, if /3 ~ /3' ~ 0.5 and if one assumed that 
the curves passed on opposite sides of the lens, then from equation (2.2), the mass 
would be Af ^ 1 X 10~^ Mq. If the mass (and Finstein ring) were this small, it 
would be easy to tell as I show below. However, since such events are rare, I do 
not consider them in detail here. Another class of rare events is when <^ 1. 

In this case it may be that |A/3| = A/3_|_. 

I focus on the typical case where |A/3| = A/3_. The main problem is dis- 
tinguishing between the ibA/3_ solutions. This degeneracy can be broken by the 
effect of the Farth's orbit. The ratio of light curves will oscillate annually as the 
Farth-source line of sight moves nearer and further from the projected MACHO 
position (see Fig. 2 in Gould 1992). The time of the largest positive fiuctuation 
differs by about 6 months between the -|-A/3_ and — A/3_ solutions. The differ- 
ence between the curves for the two solutions has an amplitude ~ 2AU//3re, i.e.. 



13 



~ 0.4% at the edge of the volume for /3 = 1. Hence they could be distinguished 
at the 3(7 level if 1% photometry were carried out once per weak. This required 
level of accuracy determines the minimum characteristics of the telescope. At any 
given time there will be ~ 30 events that must be followed. Assuming a typical 
apparent magnitude of V ~ 23, it would take ~ 60hr/week of 1 meter telescope 
time (including overhead) to get 1% photometry on these objects each week. An 
additional 60 hr/week would be required to get 3% photometry on the remaining 
1900 program objects every other month. Thus the telescope must be close to Im 
in diameter. 

The motion of the Earth can easily distinguish between the Al3± solutions in 
the case mentioned above when A/3_ ^ 1 and A/3_|_ ~ 0{1). The latter solution 
has an Einstein ring which is smaller by a large factor. Hence the amplitude of the 
annual oscillation is much larger. Since the experiment is designed to detect the 
smaller oscillations at the 3(7 level, this large effect would clearly stand out. 

6. Conclusions 

An experiment to measure the transverse velocities of galaxies with an accuracy 
of ~ 60kms~"^ at a rate of ~ 3yr appears to be generally feasible. A Im telescope 
in a Neptune-like orbit would be required. The experiment would yield information 
on large scale structure that would largely be free of systematic errors and so would 
provide a powerful test of theories. In addition, it could be used to measure the 
Hubble parameter to ~ 5%, to measure the stellar mass-to-light ratio of galaxies, 
and to probe the central regions of quasars. The largest uncertainty in designing 
the experiment is the effect of quasar variability. A variability study of quasars 
found in the MACHO observations (Alcock et al. 1993) now underway should allow 
detailed modeling of this effect. 

Acknowledgements: I would like to thank K. Ereeman, P. Osmer, B. Peter- 
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